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ft ; ABSTRACT 



IGR J08408-4503 is a supergiant fast X-ray transient discovered in 2006 with a confirmed 
association with a 08.5Ib(f) supergiant star, HD 74194. We report on the analysis of two 
outbursts caught by SwiftlBPX on 2006 October 4 and 2008 July 5, and followed up at softer 
energies with Swift/XRT. The 2008 XRT light curve shows a multiple-peaked structure with 
an initial bright flare that reached a flux of ~ 10~ 9 erg cm~ 2 s _1 (2-10 keV), followed by two 
equally bright flares within 75 ks. The spectral characteristics of the flares differ dramatically, 
with most of the difference, as derived via time-resolved spectroscopy, being due to absorbing 
column variations. We observe a gradual decrease of the Nh, derived with a fit using absorbed 
power law model, as time passes. We interpret these Nh variations as due to an ionization 
effect produced by the first flare, resulting in a significant decrease in the measured column 
| density towards the source. The durations of the flares, as well as the times of the outbursts 

suggest that the o rbital period is ~3 5 days, if the flaring activity is interpreted within the 
framework of the Sidoli et al. 1 (120071) model with the outbursts triggered by the neutron star 
passage inside an equatorial wind inclined with respect to the orbital plane. 
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1 INTRODUCTION 

IGR 108408^503 belongs to the class of the supergiant fast X- 
ray transients (SFXTs), which are transient sources in high mass 
X-ray binaries (H MXBs) associated with blue supergiant com- 
panions (see e.g., ISmith et af I l2006h . Several members of this 
class have been discovered since the launch in 2002 of the 
INTEGRAL satellite, thanks to its monitoring of the Galac- 
tic plane performed with a large field of view and a good sen- 
sitivity at hard X-rays. Indeed, several of these transients are 
highly absorbed in X -rays, and previous missions had failed to 
detect them (see e.g. ISguera et alj 120051) . Two members of this 
class are X-ray pulsars: IGR Jl 1215-595 2 dSwank et alj2007h and 
AX J1841.0-0536/IGR J18410-0535 teamba et al .11200 lh . The 
X-ray spectra of SFXTs are very similar to those typically dis- 
played by accreting X-ray pulsars, with a flat power law spectrum 
below 10 keV and a high energy cut-off around 10-30 keV (see 
iRomano et alfeooljl for the first truly simultaneous wide-band X- 
ray spectrum of a SFXT, from 0.3 to 100 keV). This spectral simi- 



larity suggests that the other members of the new class of HMXBs 
probably host a neutron star, as well. 

The emerging picture of the SFXTs as a class, especially 
thanks to m ore sensitive monito ring performed with Swift of four 
SFXTs fsee lSidoli etail J2008b) for the strategy of this on-going 
campaign], is that they sporadically undergo bright flares, up to 
peak luminosities of 10 36 -10 37 erg s _1 , with a duration of a few 
hours for each single flare. This bright flaring activity is actu- 
ally part of a muc h longer dur ation outbu r st even t, lasting days 
fsee IRomano et alj J2007I) and ISidoli etafl d2008al) l. Their long- 
term emission is characterized by a much lower intensity and flar- 
ing activity, with an average luminosity of 10 33 -10 34 erg s _1 
dSidoli et al.[2 008b). The true quiescence has been rarely observed, 
and is characterized by a very soft spectrum, with a lum inosity as 
low as 10 32 erg s" 1 din't Zanoj|2005l . lLevder et"aT]|2007h . 

IGR J08408-4503 was discovered on 2006 May 15, during a 
short bright flare with a duration of 900 s, reaching a peak flux of 
250 mCrab in the 20-40 keV energy band dGotz et al.l2006l) . Anal- 
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ysis of archival INTEGRAL observations of the source field led 

to the discovery that IGR J08408^1503 is a recurrent transient, be- 

i ■ it — ' — -\ 

ing also active on 2003 July 1 (Mereghetti et al. 2006). The refined 

position with Swift/XKT dKennea & Campanall2006h . during a Tar- 
get of Opportunity observation (ToO) performed about one week 
after the discovery flare, allowed the confirmation of the associ- 
ation with a Q8.5Ib(f) supergia nt star, HD 74194, located in the 
Vela region jMasetti et alj2006T) at a distance of about 3 kpc. Opti- 
cal spectroscopy of HD 74194 just a few days after the 2006 May 
15 flare revealed variability in the Ha profile and a radial velocity 
variation in the He I and He II absorption lines with an amplitude of 
about 35 km s - . The study of three addition al flares observed with 
INTEGRAL and Swift bv lGotz et all J2007t) led these authors to 
suggest that the orbital period of IGR J08408^503 is probably of 
the order of 1 yr, the spin period of the putative neutron star could 
be of order of hours, and the surface magnetic field is probably 
around 10 13 G. 

In this paper we report on the detailed data analysis of both the 
2006 October 4 and the 2008 July 5 outbursts caught by Swift/BKT 
and followed up at softer energies with Swift/XKT. 
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2 OBSERVATIONS AND DATA REDUCTION 

IGR J08408^1503 triggered the Swift/B AT on 2006 October 10 
14:45:43 UT (image trigger= 232309) and on 20 08 July 5 21:14:15 
UT (image trigger=3 16063. IWard et al] [2008). respectively. On 
both occasions, Swift slewed to the target, allowing the narrow 
field instruments (NFIs) to be pointing at the target ~ 1885 s and 
~ 130 s after the trigger, respectively. We note that this source 
is normally not detectable in the BAT Transient Monitor, but has 
shown a few untriggered outbursts over the past three years. Ta- 
ble Q] reports the log of the Swift/XKT observations used for this 
work. 

The XRT data were processed with standard procedures 
(XRTPIPELINE vO.11.6), filtering and screening criteria by using 
FTOOLS in the HEASOFT package (v.6.4). We considered both 
WT and P C data, and selected event grades 0-2 and 0-12, re- 
spectively jBurrows et alj|2005l) . When appropriate (observations 
00232309000 and 00316063000), we corrected for pile-up by de- 
termining the size of the PSF core affected b y comparing the ob- 
served and nominal PSF ( Vaughan et al. 2006), and excluding from 
the analysis all the events that fell within that region. Ancillary re- 
sponse files were generated with XRTMKARF, and they account for 
different extraction regions, vignetting, and PSF corrections. We 
used the spectral redistribution matrices vOlO in CALDB. 

The BAT always observed IGR J08408-4503 simultaneously 
with XRT, so survey data products, in the form of Detector Plane 
Histograms (DPH) with typical integration time of ~ 300 s, are 
available. Furthermore, event data were collected during observa- 
tion 00316063000. The BAT data were analysed using the standard 
BAT analysis software distributed within FTOOLS. BAT mask- 
weighted spectra were extracted over the time intervals simultane- 
ous with XRT data (see Sect. |3.2| l. Response matrices were gener- 
ated with BATDRMGEN using the latest spectral redistribution ma- 
trices and the default computa tion method (METHOD =MEAN) ac- 
cessing calibration parameters (Markwardt et al. 2007) that are av- 
eraged over the entire BAT field of view (FOV). For our spectral fit- 
ting (XSPEC vl 1.3.2) we applied an energy-dependent systematic 
error, amounting to 4— 5% through all the 15— 90keV energy band 
(Markwardt et al. 2007). Inspection of the BAT FOV shows that no 
source substantially brighter than IGR J08408-4503 was visibile 



Figure 1. XRT light curves of the 2006 October 4 outburst background- 
subtracted and corrected for pile-up, PSF losses, and vignetting, rebinned 
with a minimum of 30 counts per bin: a) \^\ keV band; b) 4-10 keV energy 
band, c) 4-10/1^1 hardness ratio. 
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Figure 2. BAT light curves of the 2008 July 5 outburst, rebinned at a 
S/N = 10. a) 15-25 keV band (units of counts s" 1 det" 1 ); b) 25-50 keV 
band (counts s — 1 det -1 ); c) 25-50/15-25 hardness ratio. 
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Table 1. Observation log. 



Sequence 


Instrument 


Start time (UT) 


End time (UT) 


Exposure 81 




/Mode 


(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 


(s) 


009^9^0QOOO 




7006-10-04. T 5- 1 7- T ^ 

1 IMn U . 1 ' .11) 


7006-1 0-04 1 6-^ 1 ■ 1 ? 


180 




XRT/PC 


7006-10-04 15-70-08 


7006-10-04 1705-56 


1472 


00232309001 


XRT/PC 


2006-10-04 18:27:05 


2006-10-04 18:42:58 


951 


00316063000 


XRT/PC 


2008-07-05 21:18:23 


2008-07-06 18:29:51 


29956 


00316063000 


XRT/WT 


2008-07-05 21:16:31 


2008-07-06 17:55:16 


447 


00030707003 


XRT/PC 


2008-07-08 13:26:02 


2008-07-08 18:23:18 


1978 


00030707004 


XRT/PC 


2008-07-09 18:26:26 


2008-07-09 18:43:04 


998 


00030707005 


XRT/PC 


2008-07-10 23:12:51 


2008-07-10 23:28:56 


965 


00030707006 


XRT/PC 


2008-07-11 00:55:02 


2008-07-11 02:41:57 


1308 


00030707007 


XRT/PC 


2008-07-12 20:25:43 


2008-07-12 22:14:56 


1319 


00030707008 


XRT/PC 


2008-07-14 20:28:10 


2008-07-14 23:47:51 


1203 


00030707009 


XRT/PC 


2008-07-15 04:30:10 


2008-07-15 23:51:55 


2384 


00030707010 


XRT/PC 


2008-07-16 19:03:10 


2008-07-16 23:59:57 


1474 


00030707011 


XRT/PC 


2008-07-17 22:27:20 


2008-07-17 22:38:25 


664 


00030707012 


XRT/PC 


2008-07-20 11:37:03 


2008-07-20 11:56:56 


1193 



a The exposure time is spread over several snapshots (single continuous pointings at the target) during each observation. 

during both the 2006 and 2008 observations, and we can estimate 
that possible contamination from other sources in the FOV cannot 
exceed 5 %. 

UVOT data were also collected simultaneously with the other 
instruments; however, due to the brightness of HD 74194 (LM Vel), 
the data were not usable in order to study the variability of the op- 
tical counterpart. 

All quoted uncertainties are given at 90% confidence level for 
one interesting parameter unless otherwise stated. The spectral in- 
dices are parameterized as F„ oc v~ a , where F v (erg cm -2 s _1 
Hz -1 ) is the flux density as a function of frequency v; we adopt 
r = a + 1 as the photon index, N(E) oc E~ T (ph cm~ 2 s~ x 
keV - ). Times in the light curves and the text are referred to their 
respective triggers. 



3 ANALYSIS AND RESULTS 
3.1 Light curves 

For both observations 00232309000 (2006 outburst) and 
00316063000 (2008), XRT light curves were e xtracted in 
severa l energy bands using the procedures described in lEvans et al.l 
d2007l) . Mask-tagged BAT light curves were created in the standard 

4 energy bands, 15-25, 25-50, 50-100, 100-150 keV, for obser- 
vation 00316063000, and rebinned to achieve a signal-to-noise 
ratio (S/N) of 10. Fig. Q] shows the XRT light curves for the 2006 
outburst, while Figs. [2] and [3] show the BAT and XRT light curves 
for the 2008 outburst, respectively. As shown in Fig. Q] the late 
start of the 2006 XRT observation (T+1885 s) only allowed us to 
follow the descent of the outburst. In contrast, the XRT 2008 data 
offer a much better coverage and show a remarkable second and 
third flares, which reached the same intensity as the first one. 

Figs. QJ and[3}i show the 4—10/1-4 hardness ratios. Fitting 
the 2006 XRT hardness ratio against a constant yields a value of 
0.60 ± 0.02 (xl = 1.64 for 33 degrees of freedom, d.o.f.), which 
corresponds a ~ 2.5(7 variation. In the case of the 2008 XRT data 
we see a definite trend for a softening of the spectrum as the source 
becomes brighter, which is reminiscent of what was observed, al- 
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Figure 3. XRT light curves of the 2008 July 5 outburst rebinned with a 
minimum of 30 counts per bin: a) 0.3-1 keV band, b) in the 1^4 keV band; 
c) in the 4-10 keV energy band; d) 4— 10/1 ^4 hardness ratio. 



beit with lower stat istics, during the 20 08 March 19 outburst of 
IGR J16479-4514 jRomano et afll2008h . No significant variation 
in the 2008 BAT hardness ratio is evident (Fig.ffJ;). Fitting the hard- 
ness ratio as a function of time to a constant model we obtain a 
value of 0.58 ± 0.11 (\l = 0.56 for 39 d.o.f.). 

A timing analysis was performed on WT data, after having 
converted the event arrival times to the Solar System Bary centric 
frame. We searched for coherent periodicity, but found no evidence 
for pulsations in the range 10ms-100 s. 
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Figure 4. Spectroscopy of the 2006 October 4 outburst. Top: simultaneous 
fit of BAT and XRT/WT data with an absorbed power law with a high en- 
ergy cutoff. Bottom: the residuals of the fit (in units of standard deviations). 
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Figure 5. Spectroscopy of the 2008 July 5 outburst. Top: simultaneous fit 
of BAT and XRT/WT data with an absorbed power law with a high energy 
cutoff. Bottom: the residuals of the fit (in units of standard deviations). 



Table 2. Broad-band spectral fits of the BAT+XRT data. 



Year N H T E c X 2 /d.o.f. 

(I0 22 cm" 2 ) (keV) 



2006 0.29Ioi4 O.3II023 11+1 0.903/69 
2008 6.7+. 2 ;° l-36to;56 > 14 1-191/45 



3.2 Broad-band spectra 

For the 2006 October 4 outburst we extracted one BAT spectrum 
nearly simultaneous with the XRT/WT data set and performed a 
broad-band fit to an absorbed power law with an exponential cut- 
off. The XRT data were rebinned with a minimum of 20 counts 
per energy bin to allow \ 2 fitting, while the BAT data were re- 
binned to a S/N > 5. We fit in the 0.3-10 keV and 15-80 keV en- 
ergy ranges for the XRT and BAT data, respectively, with XSPEC 
(vl 1.3.2). Similarly, we extracted a BAT spectrum simultaneous 
with the XRT/WT data (137-246 s) for the 2008 July 5 outburst. 
The best fit parameters are reported in Table [2] while the spectra 
are shown in Figs. [4] and [5] respectively. We note that there are 
no BAT event data during the second and third flares of the 2008 
light curve, since the source was already identified on-board as the 
current Automated Target and it did not reach the required count 
rate (twice the previous triggering rate plus 7 a) to trigger the BAT 
again. We also note that the presence of a cutoff is required by the 
data, since the trend in the residuals of a simple power-law model 
clearly indicates the presence of curvature in the spectrum. 

3.3 Time-resolved spectra 

Upon examination of the light curve and the available counting 
statistics, we selected different time bins over which we accumu- 
lated spectra. For the 2006 data, we extracted one spectrum in 
WT mode and two in PC mode (observations 00232309000 and 
00232309001 in TableQ). For the 2008 data, our time selections in- 
clude (see Fig.(6^,c, and Table[3): i) the initial WT data (137-246 s, 
hereafter WT1); ii) a low phase, with a mean CR ~ 3.5 x 10~ 2 



counts s _1 (4762-30180 s, orbits 2-6 of observation 00316063000, 
PCI); Hi) an intermediate phase, with a mean CR ~ 0.1 counts s -1 
(33697-64903 s, orbits 7-12, PC2); iv) the second flare, (68429- 
70690 s, orbit 13, PC3); v) the rise of the third flare (74208- 
74462 s, orbit 14, WT4); vi) the decay of the third flare (74463- 
76477 s, orbit 14, PC5). The spectra were rebinned with a minimum 
of 20 counts per energy bin to allow \ 2 fitting (see Fig.[6j;). How- 
ever, the Cash statistic JCashlll979h and spectrally unbinned data 
were used, due to the low counting statistics, to fit the PCI spec- 
trum (the unbinned spectrum is not shown Fig.0. Each spectrum 
was fit in the 0.3-10 keV energy range, adopting an absorbed power 
law model. The best fit parameters are reported in Table[3] Remark- 
ably, we obtain significantly different fitting parameters for the data 
of the three flares, despite quite similar flux levels. Indeed, for the 
first one (WT1) we obtain a photon index Twti = 1.42+q' 4 ^ and 
an absorbing column density of ./Vhwti = (6.63+. 2 'gg) x 10 22 
cm -2 , while the second and third ones (PC3, WT4/PC5) yielded 
T W T4 = 0.65+°' 15 and r PC5 = 0.78+.°;:^, ^ an absorbing col- 
umn density of AWt4 = (0.50±g;i|) x 10 22 cm" 2 A^pcs = 
(0.64lo 2 g) x 10 22 cm~ 2 . Fig.[7](bottom) shows the contour levels 
for the column density vs. the photon index, which we created to in- 
vestigate these spectral variations. As can be seen, the photon index 
range is restricted to F ~ 1 ± 0.5, while the measured absorbing 
column varies by an order of magnitude with a definite trend for 
a decreasing A r H as time goes by. It must also be noted that the 
2006 spectra, taken about 2000 s from the trigger, are quite similar 
to these 'late' (PC3 onwards) spectra. We note the presence of a 
bump in the residuals at soft energies, around 1 keV (see Fig. [5j- 
This feature, which has a variable contribution to the spectra and 
is stronger in WT1, is probably due to a combination of a slight 
spectral variability throughout the time interval during which the 
spectrum was accumulated (see Fig. [3ji and [6p), and some resid- 
ual ca libration uncertainties below 2 keV in WT data dGodet etall 
2008). We are however confident that the latter do not weaken our 
findings of a decreasing Nh, since their effect is in the opposite 
direction and we shall not discuss this feature any further. 
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Figure 6. Time selections for XRT spectroscopy of the 2008 July 5 outburst, 
(a): XRT 0.3-10 keV light curve; the intervals over which we performed 
spectroscopy are marked with vertical lines and labels. Filled circles are 
WT data, empty circles PC data, (b): 3-10/0.3-3 keV hardness ratio, (c): 
detail of (a) to show the second and third flare; 



Figure 7. XRT time-selected spectroscopy of the 2008 July 5 outburst: 
(top): Spectra extracted in the different intervals, (bottom): 68%, 90% and 
99% confidence level contours adopting an absorbed power-law model. The 
labels refer to the time intervals shown in Fig.|6]and Table|3] 



4 DISCUSSION AND CONCLUSIONS 
4.1 Multiple flaring behaviour 

In this paper we report on two outbursts of the supergiant fast X- 
ray transient IGR J08408-4503 observed by Swift on 2006 October 
4 and 2008 July 5. During the 2008 outburst, a multiple-flare be- 
haviour is observed so clearly in a SFXT, and for the first time with 
a large dynamic range, thanks to the good sampling of the source 
light curve. The 2008 outburst emission is composed of three main 
bright flares, each exceeding 10 36 erg s _1 , with a duration of a few 
hundred seconds. The flares are separated by ~70 ks (between the 
first and second flare) and ~4 ks (second to third flare). 

Structured and complex peaks during outburst flares have 
already been observed in other SFXTs [see e.g. the flaring ac- 
tivity observed duri ng the brightest ph ase of the outburst in 
IGR Jl 1215-5952 dRomano et at] 120071). or the seve ral flares 
caught in XTE J1739-302 with Swift dSidoli et al.l2008ij) . and with 
INTEGRAL at different epochs dBlav et al.l2008l) l. However it has 
never been observed spanning almost three orders of magnitude in 
flux, as found between the first two flare peaks in 2008 (WT1 and 
PC3 data set in Fig. [6] exceeding 10 counts s _1 ) and their inter- 
flare emission (PCI, which goes down to ~ 0.01 counts s~ ). The 
hardness ratios (Fig. [3jl and |6p) show an anticorrelation with the 



source intensity, and our time-selected spectroscopy demonstrates 
that this evolution is produced by a change in the absorbing col- 
umn density, instead of a change in the source spectrum (Fig.[7]l. In 
particular, there is evidence for a progressively decreasing column 
density with time (also during the low intensity emission of the part 
PC2 of the light curve). This is strongly suggestive of a ionization 
effect produced by the first flare, resulting in a significant decrease 
in the measured column density towards the source. 

During the first flare in 2008 BAT events were also available 
and a broad band spectrum could be extracted (see Table [3); how- 
ever, during the second and third flares BAT did not trigger (see 
Sect. l3.2t and the source was only detected at the ~ 6a level in sur- 
vey data, so no broad band coverage is available. The broad band 
spectra during the two flares in 2006 and 2008 are significantly 
different, especially as far as the absorbing column density is con- 
cerned, which was about 20 times higher in the first flare in 2008 
than during the 2006 flare [where the absorption is consistent with 
the interstellar value of 3 x 10 21 cm -2 derive d from the optical e x- 
tinction to the optical counterpart HD 74194 fcevderetal]|2007l) l. 
We note that the 2006 flare spectrum was extracted about 2000 s 
after the BAT trigger; Swift made a delayed slew to the source, too 
late to catch the peak of the first flare which triggered the BAT. This 
implies that the low column density observed during the 2006 flare 
is probably so low when compared with the first 2008 flare (WT1, 
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Table 3. Spectral fits of the XRT data. 



Year 


Spectrum 


Times 

(s) 


(10 22 cm" 2 ) 


r 


Flux a 


X 2 /d.o.f. 


2006 


WT 


1893-2064 


T7+ - 17 
U " i ' -0.14 


0.67tS;i? 


5.7 


0.902/42 




PC0OO 


2064-8413 


n 0,0+0.35 
u ' au -0 .28 


0.8St° f 2 


1.4 


1.477/31 




PC001 


13283-14234 


47+ ' 41 
u -*' -0.27 


49+0-30 
u -^ a -0.26 


1.5 


0.554/19 




all 


1893-14234 


0-57±g;« 


gg+0 .13 
u - oa -0.12 




1.118/96 


2008 


WT1 


137-246 


6.63±?j| 


1 49+ ' 52 


13.4 


1.255/35 




PCI 


4762-30180 


n.osl^ 2 


fi c i +0 ' 75 


0.086 


578.4 (55.78%) b 




PC2 


33697-64903 


I KA+0.40 
L -°*-0.33 


99+0-20 


0.14 


1.423/70 




PC3 


68429-70690 


"V7+ - 20 
U '°'-0.16 




1.8 


1.368/70 




WT4 


74208-74462 


50+ ' 24 




5.3 


1.371/62 




PC5 


74463-76477 


64+ ' 24 




1.3 


0.829/44 



a Unabsorbed 2-10 keV flux in units of 10 10 erg cm 2 s 1 , 

b Cash statistics Cstat, and percentage of realizations (10 4 trials) with statistic > Cstat. 



taken only 137 s after the BAT trigger), only because of the ion- 
ization effect caused by the flare that triggered the BAT (and not 
observed by XRT). 

This suggests that, during extended bright flaring activity, it is 
not possible to directly derive the mass of the clump which is sup- 
posed to produce the enhanced accretion lu minosity, in the frame- 
work of the clumpy stellar wind model liWalter & Zurita Heraj 
l2007h , because part or all of the accreting material could have been 
ionized by a previous bright flare (as in 2006, where the presence 
of a precursor flare is deduced from the fact it triggered the BAT). 

The broad-band spectrum during the 2006 flare is well fitted 
with a power law with a high energy cutoff at ~ 1 1 keV. This is 
compatible with a su rface magnetic field of a typical X-ray pulsar 
JCoburnetaT]|2002l), around 10 12 G, much lower th an that sug- 
gested bv lBozzo et alj J2008l) . or bv lGotz et all J2007I) . 

The inter-flare emission is highly variable, showing a trend 
with an average decreasing flux after the first flare (for about 10 ks), 
then continuously increasing by an order of magnitude (with a low 
intensity flaring activity, see PC2 part) up to the second flare. The 
spectrum in this lower intensity inter-flare emission is well fitted 
with a hard power law with a photon index of ~1, which implies 
that accreti on is still p r esent. The true quiescence has been ob- 
served by iLevder et alj d2007t) . and displays a much softer spec- 
trum and a lower luminosity around ~10 32 erg s _1 . This inter-flare 
emission is probably due to the accretion of inter-clump medium, 
implying a density (and/or a velocity) contrast in the supergiant 
wind of at least 10 3 , which is not unusual in line-driven winds. In- 
deed, in recent years evidence has been accumulated that hot stellar 
winds are not smo oth, as previous l y thought, but highly inhomo- 
geneous [see, e.g.. lOskinova et all J2007I) and references therein], 
thus implying important consequences, for example, in the mass 
loss rate determination, and also in the stellar evolution. 



4.2 Clumpy winds model: an estimate of the neutron star 
distance from the supergiant companion 

This new Swift data set interestingly allows to accurately determine 
the flare duration, particularly of the second (PC3 data) and third 
flares (WT4/PC5). Fitting their light curves with Gaussian profiles, 
we derive a FWHM of 370±20 s and 500±30 s, respectively. These 
durations can be used to derive information about the distance of 



the neutron star from the supergiant companion during the out- 
burst, since the clump sizes increase with the distance from the 
supergiant star, as follows. The size of a clump is determin ed by 
the balance pressu re equation. Following lLucv & White! dl980h and 
iHowk et al.1 feOOCl) . the average density of a clump is: 



a« + CpUJ 2 



(1) 



where p w (r) is the density profile of the homogeneous (inter- 
clump) wind, a w and a c are the inter-clump wind and the clump 
thermal velocity, respectively: a 2 = kTm and at = , where 
k is the Boltzmann constant, T m and T c are the temperatures of 
the inter-clump wind and of the clump, respectively. The constant 
C p — 0.29 accounts for the confining effect of the bow shock pro - 
duced by the ram pressure around the clump dLucv & Whiteri 980), 
while lu is the relative velocity between the wind and the clump 
(u) = v w — v c ). 

Wind structure simulations show that the average relative ve- 
locity ijj is ~ 5 x 10 7 cm s _1 , the average clump temperature is 
T c ~ 10 5 K, the average t emperature of the homogen eous inter- 
clump wind is T w ~ 10 7 K dRunacres & Owockil2005l) . Adopting 
these values for T c , T w , u, and /i — 1.3 (for solar abundances), we 
obtain: 



a 2 , + Cpto 2 



200 



(2) 



Since the density radial profile p w (r) of the homogeneous 
inter-clump wind is: 

Rq B v(Rob) 



p w (r) = p w (Rob)'- 



r 2 v(r) 



(3) 



where Rob is the radius of the supergiant (i?os=23.8 R©, 
IVacca et al]| 19961) , then from equations 0, {2j, ^3) we obtain: 



pc{r) = Pc(Rob) 



R 2 ob v(Rob) 



(4) 



r 2 v(r) 

where Pc(Rob) = Pw(Rob) x 200. Assuming a spherical geom- 
etry for the clumps, it is possible to obtain the expansion law of the 
clump from equation $4$: 



R c (r) = R c (R OB ) [ iP^- 
1 Rqb v o 



1/3 



(5) 
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where vo = v(Rob) is the initial velocity of the clump at the 
surface of the supe rgiant, which can be assumed vo=10 km s _1 
dCastoretalJll975h . For v(r) we consider the standard /3 law for 
the velocity profile in line-driven winds, with a terminal velocity of 
1900 km s _1 and an exponent /3=0.8. 

The flare duration t is linked with the thickness of the clump, 
therefore 



2R C 

Vc 



2R c (R OB ) 



r 2 v(r) 
r ob v o 



1/3 



(6) 



where we assumed R c (Ros) ~ 2x 10 9 cm; thus, flare durations of 
the order of ~500 s can be produced at a radial distance of ~5 Ro b 
(or 10 13 cm), which implies an orbital period of ~35 days, assum- 
ing a companion mass of 30 Mp. This estimate is in agreement with 
that obtained by llxvder et al. l 2007t). bu t is significantly shorter 
than what calculated by GQtz et alT i 20071) . 



4.3 The orbital period and the outbursts double periodicity 

We note that the first three outbursts reported in literature (on 52821 
MJD, 53870 MJD and 54012 MJD) are indeed spaced by a multi- 
ple of 35 days (30 and 34 cycles, respectively, after the first ob- 
served outburst on 52821 MJD). On the other hand, the 2008 out- 
burst (54652 MJD) does not occur strictly after an integer number 
of cycles with this assumed period. This suggested to try different 
combinations of two numbers (two sub-periods Pi and P2, such that 
P=Pi+P2=35 days), to test the possibility, suggested in Sidoli et al. 
(2007), that the SFXTs outbursts occur with a double periodicity, 
when the neutron star crosses the equatorial wind of the companion 
twice along the eccentric orbit (see their fig. 9c). In this hypothesis, 
the times of the n t h outburst is given by t„.=to+(ni Pi)+(n2 P2). 
Among the several couples of numbers we tried (with sum equal 
to 35 days), we find that Pi = l 1 days and P2=24 days well account 
for all the four reported outburst times (calculated with m=n2=30, 
m=n2=34 and with ni=53 and n2=52, respectively for the three 
outbursts, assuming to=52821 MJD), within an uncertainty of 1 
day. 



4.4 IGR J08408-4503 in the SFXT context 

The multiple flares we observe here, reaching similar peak lumi- 
nosities, cover about one day of bright outburst emission and are 
then followed by lower X-ray emission in the subsequent days, 
similar to what we already observed in four other SFXTs with Swift 
(see Fig. [8). The three flares observed here so close in time, after 
a long period during which no bright flares have been observed, 
are suggestive of the fact that the neutron star is probably cross- 
ing a region of higher average density of the wind, compatible 
with the model we proposed for IGR Jl 1215— 5952 ( Sid oli et"al] 
l2007h . In this hypothesis, the outbursts in SFXTs are produced 
when the neutron star crosses an equatorially enhanced region of 
the supergiant companion, which is inclined with respect to the or- 
bital plane. The different flares then are probably produced by in- 
homogeneities in the structure of the wind, which, in our proposed 
model, is anisotropic, showing a preferential plane along the equa- 
tor of the supergiant companion. This could explain why the flares 
happen in a preferential region of the neutron star orbit, display- 
ing a particularly high X-ray flare rate. New observations during 
the next outbursts are crucial in order to test this hypothesis and 
to confirm whether multiple outbursts are indeed a usual property 
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Figure 8. Light curves of the outbursts of SFXTs followed by Swift/XKT 
referred to their respective triggers. We show the 2005 outburst of 
IGR J16479-4514 jSidoli et al.ll2008bl). w hich is more complete than the 
one observed in 2008 ^Romano et alj .2008). The IGR Jl 1215-5952 light 
curve has an arbitrary start time, since the source did not trig ger the BAT 
(the observations were obtained as a ToO; Romano et al. 20 03). The data 
on XTE J1739-302 and IGR J17544-2619 are presented in ISidoli et alj 
l 2008a). Note that where no data are plotted, no data were collected. For 
clarity, we drew dashed vertical lines to mark each day (up to one week) 
since the trigger. 



of this (and other) SFXTs, and if there is indeed an orbital phase 
dependence of the X-ray flare rate. 
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